Anisotropic Elastic Model for Short DNA Loops 
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Effect of bending anisotropy on a planar DNA loop, using energy minimization and neglecting 
entropic effects, is studied. We show that the anisotropy results in polygonal shape of the loop and 
increasing the anisotropy makes the edges sharper. Calculating the energy of such a loop lets us to 
find effective persistence length as the geometrical mean of hard and soft rigidities, which is quite 
different from harmonic mean for an unconstrained long DNA. 



I. INTRODUCTION 

A long time has been passed since a straight double he- 
lical structure had been proposed by Watson and Crick 
as relaxed configuration of a DNA, but today, we know 
it can almost never be found in this form in nature. The 
macromolecule is able to play its important and essen- 
tial role in the life, when it is interacting with differ- 
ent proteins which force it to bend, twist, melt and/or 
pack 1,2}. 

To understand the macromolecule's functionality, 
DNA has been examined in different length scales. Sin- 
gle molecule force-extension studies on 10 — 100// m 
DNAs H, cyclization experiments on 100 — 1000 nm 
DNAs [4] and recent experiments on sharply bent DNAs 
in loops with the lengths less than 50 nm, persistence 
length of DNA [|, 1, 0, Q . All these experiments are of 
great biological interest. Depending on relevant length 
scale of the problem, experimental results are usually an- 
alyzed in two types of theoretical models: variants of the 
continuum elastic worm-like chain (WLC) and base-pair 
steps. 

The worm-like chain models are based on the physics of 
elastic rod where the energy has a quadratic form in de- 
formations. There are two approaches to this model. The 
first one is based on Kirchhoff-like equations of balance of 
forces and torques for every segment of the rod [jjl llbTl and 
the other one is based on energy minimization (n[ • 
In these models statistical mechanics of the rod at non- 
zero temperature is also encountered 13]. A force con- 
stant or equivalently a deformation module is defined for 
bending, twisting, stretching and their coupling terms. 
In the classical WLC model, all cross terms as well as 
stretching energy are neglected and the rod is isotropic 
in bending in different directions [lij]. So the model is 
identified with just two persistence lengths assigned to 
each of bending or twisting degrees of freedom. Modi- 
fications to this model have been made by encountering 
twist-stretch coup ling, tw ist-bend coupling and stretch- 
bend coupling [U,[II[IIE[l2i[II| and estimating the free 
parameters by fitting the theory to force-extension exper- 
imental data. WLC model has also been used to evaluate 
loop formation probability, J- factor [l2|, [Hj] ■ 



Looking at microscopic structure of DNA macro- 
molecule suggests that bending toward the groove is 
easier than bending toward the backbone. Though, 
there is no direct experiment to measure the anisotropy, 
there are some theoretical postulates for double stranded 
DNA 0, H(| HH • Monte Carlo simulations have also 
confirmed that B-DNA bends more easily in the groove 
direction (roll) than in the backbone direction (tilt) [22| . 
It has been encountered in some models for describing 
force-extension and DNA cyclization experiments and 
found to have no significant improvement to the re- 
sults [ll|, |l2j, |23j , though existence of bending anisotropy 
affects the existence of twist-stretch coupling [l7| ■ Balaeff 
et al. have claimed that bending anisotropy can reduce 
the energy of loops in length scale of 76 bp by a factor of 
one third [24l |. In a recent experiment an asymmetry in 
the periodic behavior of free energy of loop formation as 
a function of loop length for 60 — 100 bp loops has been 
detected 0j- Fourier analysis of free energy shows two 
main frequencies. One with a period of ~ 10.5 bp due to 
helical shape of the double strand and one with a period 
of ~ 5.6 bp that might be a result of bending anisotropy. 

Sequence dependence and anisotropy of bending per- 
sistence length has been widely noticed in base-pair steps 
approaches, in which relative rotation and displacement 
of every two segments are defined through six parameters 
slide, shift, rise, tilt, roll and twist [25j. These sequence 
dependent parameters for individual base-pair steps have 
been determined from their standard deviation in crystal 
complexes |26j. Theoretical work has also been done on 
extracting these parameters from atomic level parameters 
via the analysis of Molecular Dynamics (MD) trajecto- 
ries (27| ■ Different simulation studies give estimations for 
roll and tilt values. Munteanu et al. have shown that for 
3 — 11 bp DNAs bending rigidity oscillates with bending 
direction and the values of roll is ^ 8 — 10 times greater 
than the values of tilt [28[ . Olson et al. estimate the ra- 
tio of the hard bending rigidity, Ai, to the soft bending 
rigidity A2 to be between 4 and 16 [25|] and between 1 
and 5 in a more recent study [2{|. MD simulations of 
17 bp dsDNAs gives A\ almost twice as A2 [30(. In a 
stack of plates simulations, Mergell et al. state that spa- 
tial constraints make roll twice as favorable as tilt 13111 . 
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Effect of sequence dependency on persistence length and 
also intrinsic curvature on loop formation problem have 
been studied analytically by Popov et al. which result in 
a wide distribution of cyclization probabilities 32 1. 

An effective persistence length, A, can be assigned to a 
rod with hard and soft anisotropic bending rigidities A\ 
and A2 (corresponding to tilt and roll) . In an analytical 
stack of plates study, done by Olson et al., a complicated 
dependence of A on detailed couplings, anisotropic con- 
stants and sequence alphabet has been shown for a free 
DNA [H. 

For a free, long and highly twisted DNA, the effective 
persistence length exactly equals the harmonic mean of 
soft and hard rigidities, 



in 



This can be easily deduced from the equipartition prin- 
cipal. The equipartition principle states that the total 
energy of a rod is equal to ^ksT times number of its de- 
grees of freedom. Bending of a segment of a rod around a 
principal axis with length equal to its persistence length 
can be considered as one degree of freedom. Therefore 
number of degrees of freedom in hard and soft directions 
is simply counted as and and the total energy 
is found to be \k B T{-^ + jfc). To find the efficient 
persistence length, the anisotropic rod is considered as 
an isotropic one. Since the isotropic rod is able to be 
bent around two similar directions, its energy is found 
to be 2(ifceT-j) in which A is the efficient persistence 
length. Having the two estimations equal to each other 
results equation [TJ This is more accurately derived by 
Kehrbaum using an averaging theory for a non-isotropic 
elastic rod with high intrinsic twist [33j . Adding geomet- 
rical constraints on a free DNA might affect the above 
relation. 

Here we are going to give an expression for effective 
persistence length of in-plane DNA loops, by energy 
minimization, when the entropic effects are neglected. 
Although an isotropic WLC loop with length smaller 
enough than its persistence length is planar, it is not 
obvious for the case of anisotropic DNA. In fact experi- 
ments [34[ and MD simulations [331 show that DNA mini- 
circles are not completely, but are almost planar. 

The rest of paper is organized as follows: Section 2 
describes the model that is used in studying the short 
DNA loops, followed by the presentation of the results in 
Section 3. Finally, Section 4 concludes the paper, while 
the overestimation of deviation of DNA local twist from 
its mean value appears in Appendix. 



II. THE MODEL 

The anisotropic elastic model of DNA represents the 
macromolecule as an elastic rod of length L, parameter- 
ized by arclength s. As the double-strand is ribbon-like, 



it is anisotropic in bending around two different direc- 
tions. On the other hand, DNA is twisted, so while bend- 
ing around a fixed axis, it should bend around each of the 
two directions, periodically. It is energetically preferable 
for the DNA to bend more around the "soft" axis and 
less around the "hard" one. To reduce elastic energy, the 
DNA also may modify its twist to have more soft bending 
along its planar path. 

For a planar DNA loop, the tangent unit vector t and 
the twist angle tp at each point contain enough infor- 
mation to parameterize DNA conformation. Actually 
an anisotropic bent DNA has a higher tendency to go 
out of the bending plane which depends on the value of 
anisotropy (Ai~ A 2 ) / {Ai+ A 2 ) . For short loops, this ten- 
dency decreases as the anisotropy becomes smaller [36j . 

Elastic energy of anisotropic rod is 
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A(s)\t\ 2 + C(tP-u;o) 2 ds, (2) 



where A(s) is the local bending rigidity, C is the twist 
rigidity, ojq is the spontaneous twist of the helix, and the 
dots indicate derivatives with respect to s. The first term 
is the usual elastic term, with s dependent persistence 
length A(s) and the second term is the energy needed to 
over- (under-) twist the DNA, which may be implied by 
the boundary conditions on ip or the anisotropic effects. 
It should be noted that in this model s dependence of 
A is not because of sequence dependence, but it is due 
to rotation of soft and hard axes. As we will show later, 
this is a source of twist-bend coupling, although there is 
no explicit coupling term in the Hamiltonian. 

DNA bending could be decomposed into two principal 
axes ei and e.2, attached to the DNA. The hard one, e\ is 
perpendicular to double strand's local plane, and the soft 
one, £2, is defined to lie in the local plane of the double 
strand, and to be perpendicular to both strands. Due 
to the helical structure of DNA, e\ and £2 rotate with 
the helix. Since t is perpendicular to t axis and lies in 
the ei — £2 plane, by its decomposition in our coordinate 
system we obtain 



E 



1 



-k R T 



(Ai sin 2 ip + A 2 cos 2 ijj) 9 2 
+ C(ii~Lo ) 2 ]ds, (3) 



where A\ and A2 are constant bending rigidities about 
the rotating axes, e\ and e~2, and 9 is the angle between 
t and a fixed, arbitrary direction in loop's plane. Us- 
ing Euler-Lagrange equation and applying corresponding 
boundary conditions, we are able to find 9{s) and ip(s) as 
well as DNA's shape and bending energy of the loop. The 
integral form of closed loop or "ends meeting" condition 



/ sin#(s)ds= / cos0(s)ds = O. 
Jo Jo 



(4) 
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Also because the DNA strands are antiparallel they can 
not switch in the ends and they should bind in phase, 
then we have ip(L) — ip(0) = 2kir with k being an integer 
that is the number of full turns of the helix. On the 
other hand to avoid any singularity in tangent vector 
changes we consider a constraint on 9 values in the ends 
by 9(L) -9(0)= 2n (no self crossing). 



III. RESULTS 

Integrating Euler-Lagrange Equations of results fol- 
lowing equations of motions: 



(A 1 +A 2 ) - (Ax -A^co^iP' 



(5) 



272 1 

" C (A 1 +A 2 ) - (A 1 -A 2 )cos2^ 



+ P, (6) 



where integral constants of 7 and should be determined 
from the boundary conditions. As it is seen in above 
equations, ip depends on s, as on the DNA's local bend- 
ing in a complicated form. This twist-bend coupling is a 
direct result of anisotropy of the model and vanishes in 
the case of isotropic rods. Even if A\ and A 2 differ by 
one order of magnitude, which is the case of our double 
stranded DNA, the coupling is very weak and s depen- 
dence of ip is negligible. Indeed, numerical studies show 
that the relative variations of ip is less than one percent, 
even if two bending rigidities differ by two orders of mag- 
nitude 36] . This leads us to consider homogeneous twist 
along the DNA and set ip(s) = los, in which 10 = 2kir/L 
is found by applying the "ends in phase" boundary con- 
dition. This is more discussed in Appendix. The above 
approximation decouples torsional part of energy from 
its bending part. 

Integrating ((5|) with considering above approximation 
gives 



9(s) = 9(0) 



7 



>VAi~A 2 



tan 



(7) 




tan(ws) + 7r 



LOS 
7T 



where the bracket means "integer part". This term is 
added to get rid of discontinuity in the tan -1 function. 
Without lack of generality, we set 9(0) = and therefore 
the condition on total bending results 9(L) = 2tt. Apply- 
ing this condition simply lets us to fix 7 = (2tt / L)\/ A\A 2 . 
These would yield the functional form of 9(s) and loop's 
shape as 




A x . 

— — tan(ws) + 7r 
A 2 



LOS 
7T 



(8) 

The above solution automatically satisfies the "ends 
meeting" condition k in above equation counts the 



number of turns of helix along the loop. To have min- 
imum torsional energy it should be fixed to the closest 
integer value to luqL/2t{, so k = round(iOoL /2ir). Thus 
the DNA is undertwisted in case Aip = ioqL — 2kn < 
and is overtwisted otherwise. 

In the case Ax » A 2 for a closed DNA loop, bending 
is not homogeneous and it is localized in "soft" parts. 
Here the loop looks more like a polygon rather than a 
circle. This is the direct effect of bending anisotropy. 
The shape of looped DNA for anisotropic model is given 
in figure [T] Increasing the ratio of Ax/A 2 , the polygonal 
shape is more visible (e.g. for Ax/A 2 > 50, the DNA 
loop will have sharp edges at the soft points). In a full 
helix turn the DNA meets the soft axis of rotation two 
times when it bends in plane, thus the number of polygon 
edges equals 2k. As the number of edges (~ helical turns) 
increases by length of the loop, the polygonal shape of 
the loop is less visible for larger DNAs. 

To find the energy of the loop, we read (Ax sin 2 ip + 
A 2 cos 2 ip) 6 = 7 from © and substitute it in 



E = -k B T 
2 B 



7 + C(iP - Lo y ds 



(9) 



As J 9ds = 2ir and ip and 7 are constant, the elastic 
energy is 



27T 2 

L 



, kvtT 

E=—k B T^fMA 2 + ^-CAip 2 . (10) 



The second term is the twist energy which is due to the 
over- (under-) twist implied by the "ends in phase" condi- 
tion on DNA loop. Because k is a step function of L, this 
term leads to a well known oscillatory behavior which is 
damping by an L _1 factor. In analogy with the bend- 
ing energy of a circular loop, (2ir 2 kBTA/L), we can read 
the effective persistence length of in-plane small loops as 
\J A\A 2 . This is different from the effective persistence 
length of anisotropic DNA in larger scales (equation |(TJ)). 



IV. CONCLUSION 

We studied the effect of bending anisotropy on planar 
DNA loops using energy minimization and neglecting en- 
tropic effects under constraints of parallel and in phase 
ends. Bending anisotropy causes anisotropy in curvature 
and twist. However, anisotropy induced a twist-bend 
coupling to the model even in the lack of explicit ap- 
pearance of such coupling in the Hamiltonian, though, 
it is small enough to be neglected in the calculations. 
The bending anisotropy results in polygonal shape of the 
loop. Increasing the anisotropy makes the edges sharper 
where the number of helical turns and hence the number 
of edges of the polygon grows with loop length. 

Energy of such a loop includes an oscillating term 
for twist energy and a bending term similar to that 
of anisotropic loop with an effective persistence length 
A = \J AxA 2 , which is different from harmonic mean of 
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L=100 bp 





FIG. 1: Above, Loop shapes for different DNA lengths. An 
exaggerated ratio of 4 1 = 100 for all illustrations has been 
used. Below, Effect o/anisotropy on sharpness of polygonal 
edges. A loop length of L — 9.6 nm for illustrations has been 
used. In all figures u> = 1.8 nm -1 (see text). 



A\ and A 2 (equation [I} for an unconstrained long DNA. 
As the geometric mean of A\ and A 2 for A\ ^ A2 is al- 
ways larger than their harmonic mean, the results shows 



that considering anisotropy in DNA bending rigidities is 
not able to explain large loop formation probability of 
DNA minicircles. Thus it seems other efforts, as like as 
bubble formation (3?| , effect of sequence dependence [13] , 
or generalized semiflexible model 13S| are more successful 
in this direction. 



APPENDIX 

We are going to give a simple numeric estimation for 
deviation of tp from its mean value. Since cos 2ij) varies 
between — 1 and 1 the maximum value of the ip dependent 
term (first term) of equation ([5]) is 



max(x) 



"7 



1 



C A, 



-4w 2 A 1 A 2 
CA 9 L 2 ' 



(11) 



Applying the typical values for a DNA minicircle (i.e. 
to = 1.8 nm -1 , C = 75 nm, L = 94 bp and bp = 
0.34 nm) following by rough estimations, A\ = 100 nm 
and A2 — 33 nm (i.e Ai/A 2 = 3 with a harmonic mean 
of 50 nm), result a max(i) ps 0.016 to 2 . So the deviation 
of ip is small comparing to the mean value and therefore 



equation 



OJ 



Now finding (max — min) of both sides of 



we have: 



7 2 1 1 4tt 2 



and as ip ; 



Alp _ Att 2 (A 1 - A 2 ) 
~ ~ L 2 Cu) 2 



1%. 



(12) 



(13) 



It should be noted that above calculation is an overesti- 
mation of the range of ip changes. As it was mentioned 
in the text, numerical study shows less than one percent 
fluctuation in the value even with considering two rigidi- 
ties different by two orders of magnitude. That means, 
considering to as a constant (-0 = ws), is not a rude ap- 
proximation. 
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